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Abstract: A common feature of retirement income products is that their payouts depend on the lifetime of 
policyholders. A typical example is a life annuity policy which promises to provide benefits regularly as long 
as the retiree is alive. Consequently, insurers have to rely on “best estimate” life tables, which consist of 
age-specific mortality rates, in order to price these kind of products properly. Recently there is a growing 
concern about the accuracy of the estimation of mortality rates since it has been historically observed that life 
expectancy is often underestimated in the past (so-called longevity risk), thus resulting in longer benefit pay¬ 
ments than ins urers have originally an ticipated. To take into account the stochastic nature of the evolution of 
mortality rates, Lee and Cartel ( 1992 ) proposed a stochastic mortality model which primarily aims to forecast 
age-specific mortality rates more accurately. 


The original approach to estimating the Lee-Carter model is via a singular value decomposition, which falls 
into the least squares framework. Researchers then point out that the Lee-Carter model can be treated as a 
state-space model. As a result several well-established state-space modeling techniques can be applied to 
not just perform estimation of the model, but to also perform forecasting as well as smoothing. Research in 
this area is still not yet fully explored in the actuarial literature, however. Existing relevant literature focuses 
mainly on mortality forecasting or pricing of longevity derivatives, while the full implications and methods of 
using the state-space representation of the Lee-Carter model in pricing retirement income products is yet to be 
examined. 


The main contribution of this article is twofold. First, we provide a rigorous and detailed derivation of the 
posterior distributions of the parameters and the latent process of the Lee-Carter model via Gibbs sampling. 
Our assumption for priors is slightly more general than the current literature in this area. Moreover, we 
suggest a new form of identification constraint not yet utilised in the actuarial literature that proves to be a 
more convenient approach for estimating the model under the state-space framework. Second, by exploiting 
the posterior distribution of the latent process and parameters, we examine the pricing range of annuities, 
taking into account the stochastic nature of the dynamics of the mortality rates. In this way we aim to capture 
the impact of longevity risk on the pricing of annuities. 

The outcome of our study demonstrates that an annuity price can be more than 4% under-valued when different 
assumptions are made on determining the survival curve constructed from the distribution of the forecasted 
mortality rates. Given that a typical annuity portfolio consists of a large number of policies with maturities 
which span decades, we conclude that the impact of longevity risk on the accurate pricing of annuities is a 
significant issue to be further researched. In addition, we find that mis-pricing is increasingly more pronounced 
for older ages as well as for annuity policies having a longer maturity. 


Keywords: Mortality modeling, longevity risk, Bayesian inference, Gibbs sampling, state-space models, life 
annuities 
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1 INTRODUCTION 

The pricing of retirement income products depends crucially on the accuracy of the predicted death or survival 
probabilities. It is no w widely documented that survival p robability is consistently underestimated especially in 
the last few dec a des dlnternational Monetary Fund ( 2012h '). To capture the stochastic nature of mortality trends, 
Lee and Carter ( 1992h proposed a stochastic mortality model to forecast the trend of age-specific mortality 
rates. 

There exists a body of literature on how to estimate the Lee-Carter model. The original approach in 
Lee and Carterl (119921) is via singular value decomp osition. To overcome the unrealistic feature of homogene¬ 
ity in the additive error term, iBrouhns et al.l (l2002h recast the model as a Poisson regression model assuming 
Poisson random variation for the number of deaths. Estimation of the model in the Poisson regression set¬ 
ting under the Bayesian framework is carried out in ICzado et al.l (120051) . Also there is a recently developed 
framework for modeling death counts with common risk f actors via credit r isk plus methodology and resultant 
estimation of the model via Monte Carlo Markov Chain in Hirz et al. ( 2015h . In this paper we focus principally 
on the class of what has become known as the Lee-Cart er models , in thi s regard another approach to estimating 
the Lee-Carter model is via state-space representation. Pedroza ( 20061) shows that the pr edictive intervals for 
foreca sting are materially wider than using the singular value decomposition method. iKogure and Kurachi 
( 2010l) adopt the state-space modeling approach and apply it for the pricing of longevity bonds and swaps. 


In this paper we aim to explore further the Bayesian state-space modeling approach and examine its implication 
for annuity pricing. Specifically, we provide a rigorous and detailed derivation of the posterior distributions of 
the static parameters and the latent process of the Lee-Carter model via Gibbs samp ling. Our assumptions on 
the pri ors on the Lee-Carter model parameters are more general than Pedrozal ( 2006 ) and Kogure and Kurachi 
( 2010l) . Moreover, a new form of identification constraint not yet recognised in the actuarial literature is 
proposed which proves to be more convenient for estimating the model using an MCMC method under the 
state-space formulation. Using the predictive distributions of age-specific death rates, we examine the impact 
of longevity risk on the pricing of annuities and demonstrate that this long-term risk is indeed a significant 
factor when accurate pricing is required. 

In Section 12] the state-space Lee-Carter model is presented together with some definitions and notation. Sec¬ 
tion |3] describes the Gibbs sampling approach to estimate the state-space Lee-Carter model. Posterior distri¬ 
butions of the static parameters and the latent process are derived in detail. Section|4|examines the impact of 
longevity risk on annuity pricing. Section|5]concludes with some remarks. 

2 LEE-CARTER MODEL 


2.1 Definitions and Notation 


In this section we briefly recall some important definitions from actuarial literature on mortality modelling 
that are required to set u p the Lee-Cyter m odel below and the pricing analysis in Section |4| We follow 
Dickson et al. ( 20091) and IPitacco et al. ( 20091) . Let be a random variable representing the remaining life¬ 
time of a person aged x. The cumulative distribution function and survival function of are written as 
t(1x = P{Tx < t) and rPx = PiT^ > t) respectively. For a person aged x, the force of mortality at age x + t 
is defined as 


Px+r ■■= lim \P{T, 

h—^r\ h 




and hence rPx = exp (— fj Px+s ds). Note that the survival probability function has the following important 
property: r+uPx = tPx+u x uPx- Let fx{t) be the density function of T^, then from ((Tj) we see that xQx = 
Jo fxis) ds = /p sPx Px+s ds. The central death rate for a a:-year-old, where a: € N, is defined as 

ife fn sPxPx+sds 

TUx '■= — - = —I - (2) 

Jo sPx ds /p sPx ds 


which is a weighted-average of the force of mortality. Under the so-called piecewise constant force of mortality 
assumption, that is px+s = Px where 0 < s < 1 and a: S N, we have, from (|2]), rux = Px and hence 
iPx = Moreover, the maximum likelihood estimate of the force of mortality p^ (and hence Ax) is 

given by p^ = P>xlPx = rhj; where Dj, is the number of deaths recorded at age x last birthday and the 
exposure-to-risk Ex. is the total time lived by people aged x last birthday, during the observation year. Note 
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that Ex is often approximated by an estimate of the population aged x last birthday in the middle of the 
observation year. 

2.2 The Lee-Carter State Space Model 


Based on the definitions described above, we now discuss the work of iLee and Carted (Il992h who pro¬ 
posed a stochastic mortality model specifically for forecasting age-specific central death rates ruxt, where 
X = xit .. ,Xp and t = 1,...,n represent age and year (time) respectively. The model assumes that the log 
central death rate, yxt = In tnxt, is governed by the following equation 


y^ = a +f 3 Kt +Et, et ~ N(0,cr^lp) 

where = {yx^t, ■ • ■, yx^t)', a = {ax ^,. • ■, 

by p identity matrix and N(.,.) denotes the Gaussian distribution. iLee and Carted (Il992h estimate the model 
Q via singular value decomposition and subsequently assume that the unobserved latent time trend denoted 
by Kt satisfies the following linear dynamics 


( 3 ) 

) . /3 = if^xi ^ ; Pxx ) j ~ (^xi 11 ■ ■ •! ^Xpt) . Ip is the p 


Kt — Kt-i + 0 + LOt, Wt ^ N(0, tr^) 


( 4 ) 


where Et and ujt are independent. The parameters 0, ctJ are then estimated using standard econometric tech¬ 
niques. In this form the Lee-Carter model is, however, not identifiable since the model Q is invariant up to 
some linear transformations of the parameters: 


yt = a. + (3Kt + Et = OL 


'/3c+ ^ [{nt - c)d) 
a 


■ £( — Q: + ^Kt -f Et 


( 5 ) 


where a = a -f /3c, /9 = ^ and kt = {Kt — c)d. To overcome this identification issue. iLee and Carteii(ll992h 
introduced the following constraints 


^ /3a: = 1, ^ Kt = 0 




(6) 


and c = X]r=i 


to ensure that the model becomes identifiable since, by setting d = J2x=x - 
T,7=xt Px = 1 and Y,t=i = 0- 

IPedrozal ( 2006 ) suggests that we can in fact combine the processes and Kt into one dynamical system 


yt = a +fdKt + Et, Kt = Kt-i + 0 + uJt, where £t ~ N(0, erf Ip), Wt ~N(0, erf) 


( 7 ) 


resulting in a state-space representation of the Lee-Carter model and estimate ttt and model parameters jointly. 

Note on Estimation via Lee-Carter Framework: although Lee-Carter model is expressed in a state-space 
formulation, given this form of the identification constraints it is not readily amenable to standard state-space 
estimat ion procedures s ince t he constraint is expressed on the path-space of the latent process. Consequently, 
this led Lee and Carterl l 1992 ) to develop an alternative estimation procedure where the first step in the estima¬ 
tion does not depend on the dynamics of Kjand utilises a singular value decomposition (SVD) approach. Then 
the evolution of Kt is specified after the SVD procedure is performed, meaning that although the state space 
structure is specified, this form of the representation is not exploited in the estimation of the model trend Kt or 
the model parameters. We will demonstrate in this paper how to change the identification constraints so that 
standard filtering based state-space model estimation procedures can be utilised. 


2.3 Lee-Carter model in ARIMA Time Series Form 


We also note that, at least when one doesn’t consider the identification constraints, the Lee-Carter model is 
a simple linear dynamic model. Hence, we also highlight that this model can be rewritten in the form of an 
ARIMA structure via a Local Level formulation where we denote rjt •= cx. + f3Kt and ht := 6 + Wt- One 
can then rewrite the state-space form where each age x is an ARIMA(0,1,1) structure as Zxt ■= V/Fxt = 
Vrjxt + ^^xt with a simple closed form expression for the auto-correlation function given by 


PZp (k) 


l{k) 

7(0) 


’ k — 1, 

0, fc > 1. 


( 8 ) 
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Suggesting that one can also perform estimation on the unconstrained form of the model via estimation based 
on the autocorrelation, though these would need to be modified subject to identification constraints. This 
would again complicate the estimation, suggesting the need to try to find alternative identification constraints 
that are more applicable to these standard estimation approaches. 


3 BAYESIAN INFERENCE FOR LEE-CARTER MODEL IN STATE-SPACE FORM 


Pedrozal ( 2006h and Kogure and Kurachi ( 201(lll both consider Bayesian formulations of the Lee-Carter model 


which allows the joint estimation of Kt and model parameters. However, under their formulation they again 
work with the identification constraints of ® which are not obvious to use when designing efficient Monte 
Carlo procedures such as an Markov Chain Monte Carlo (MCMC) procedure. Such identification constraints 
will lead to difficulties in designing the proposal of the MCMC and difficulties in achieving suitable acceptance 
rates for the resultant Markov chain, resulting in high variance in estimates of mortality rates. 


Additionally, although these authors work in the Bayesian setting, their derivations of the posterior distribu¬ 
tions are not fully described. In the following we derive the posterior distributions of the parameters and the 
state process of the Lee-Carter model under our extended Bayesian framework. 


3.1 Lee-Carter model: New Identification Constraints and Bayesian Formulations 


In this article, we suggest an alternative new formulation of the identification constraints required which we 
believe is simpler and more readily applicable to most Monte Carlo based procedures such as MCMC and 
filtering methods such as Kalman Filter and Sequential Monte Carlo. This has the key advantage that for a 
given computational effort we can design efficient MCMC samplers with lower variance and therefore result 
in more reliable estimates of mortality rate. Our formulation of the identification constraints are given by 
simply setting = constant, and /3xi = constant. Such a choice is a valid identification constraint since if 
one of the elements of each a and {3 are known, then a non-trivial linear transformation in Q is not allowed 
anymore; that is, we must have c = 0 and d = 1. 


Under the Bayesian approach, we aim to obtain the posterior density 7r(Ko:n, of statesQ 

well as the parameters, := o'ei given the observations 


Ko:„ as 


Note that Ux^^ and , 


'Xl 


are assumed to be known constraints. Under such a Bayesian formulation, it is standard to utilise a MCMC 
pr ocedure to sample from 7r(Ko:n, discussions on such procedures in risk and insurance settings 


Cruz et aP ( 2015 ). 


In this paper we explain an efficient and suitable sampling framework for actuarial applications which utilises 
the state-space Lee-Carter structure, in particular the fact that it is a linear Gaussian model, as well as the new 
constraint formulation we introduce. Under this model we develop an efficient approach involving a combined 
Gibbs sampling conjugate model sampler for the marginal target distributions of the static model parameters 
along with a forward backward Kalman filter sampler for the latent process 

A sample of the targeted density is obtained via Gibbs sampling in two steps: (1) Initialise ^ (2) For 

i = 1,... ,N, first draw fromthen draw from7r(tl'|«;g*^, Vi-n)- 


3.2 Sampling from the full conditional density 7r(/to:„ |’S', „) 


Samples from the full conditional densi ty Trf/tn:,;. 1^. y^ can be obtained via the so-called forward-filtering- 
backward sampling (FFBS) procedure (iCarter and KohnIdlQQ^B . We can write 


r(K0:n|4',yi:„) = 7r(/tt ’S' , y 


( 9 ) 






where the last term in the product, 7r(Kn|’®^, Vi-n)^ distributed as N(mn, Cn) in Kalman filtering. We use 
the following notation 


Kt-i|yi:t_i ( 10 ) 
^ N(at,i?t), where at = mt-i + 0, Rt = Ct-i + al (11) 

yt|yi:t-i ^ where f t = a + f3at,Qt = (30Rt + a'^lp (12) 

Kt\yi,t ^ ^{mt,Ct), where mt = at + Rt(3'Qt^ [yt - ft),Ct = Rt - Rtf3'Qt^/3Rt (13) 

^Here ai-.t means ai,..., at. 
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to denote the distributions involved in Kalman filtering. Once we draw a sample from N(m„, C„), then 
Eq. (|9]l suggests that we can draw recu rsively and backw ardly Kt from where t = n — 

1, n — 2,..., 1, 0. It can be shown that ( Petris et ^ ( 20091) 1 

where ht = mt + CtRt+i{Kt+i-at+i),Ht = Ct-CtRt+iCt. (14) 


In summary, the EEBS algorithm consists of three steps: (1) Run Kalman filter to obtain m„ and Cn\ (2) Draw 
Kn from N(m„, Cn) and (3) Eor f = n — 1,..., 0, draw Kt from N{ht,Ht). 


3.3 Sampling from the full conditional desnity 7r(’®'|Ko:„, Vi-n) 


Sampling from the full conditional density 7r(4'|Ko:n, yi „) can be achieved by applying Gibbs sampling. 
The prior for {a^, Px,0,a‘^,al) are given by Px ~ tr^ ~ IG(de,&e), 9 ^ 

N{jlg, (f|), (T^ ^ IG(d(^, buj) where x G {x 2 , ■ ■ ■, Xp} and IG(.,.) denote the inverse-gamma distribution. It 
is assumed that the priors for all parameters are independent. In this case the posterior densities of parameters 
are of the same type as the prior densities, a so-called conjugate prior. In the following we derive the posterior 
distribution for each parameter (for ease of notation it is assumed that y = k = Kg-n, family means 
“’®' without the parameter A”): 


• Eor ax where x G {x 2 ^ ■. ■, Xp}, we have 

n 

'x{ax\y, K, oc TT{y\K, ■4')7r(K|\P)7r(aa^|’S'_Q^) oc n T^iVxtlKt, ax, Px, (rl)'K{ax) 


1 ({din + al)al - + dl - PxKt))c 


oc exp 

Hence the posterior conditional distribution of ax is given 
• Eor Px where x G {x 2 , ■ ■ ■, Xp}, we have 


byN 


’ cr2n+cr| 


5-2 n+crf 


Tr{Px\y,K,^_pJ oc Tr{y\^)Tr{K\<P)TT{Px\^-pJ o^W_'K{yxt\Kt,ax, Px,(rl)'x{Px) 

t^l 

1 ( (^1 Et + ^l)Pl - 2 (A/ 30 -e + cr| “ OLx)Kt^ Px 


« exp ^ -- 




Hence the posterior conditional distribution of Px is given by N f ^ ) 

• Eor 9, we have 

n 

TT{9\y,K,'l>^g) oc 7r(y|/«, ^)7r(/«|^)7r(6>|^^e) oc 7r(«:t|Kt_i, 0, cr^)7r(6>) 


(X exp < — 


1 / (o-|n -f al)9‘^ - 2 (/iecr^ ^2 J2^{Kt - Kt-i)) 9'' 


f2n-2 


Hence the posterior conditional distribution of 9 is given 
• Eor al, we have 


byN ( 






ain-\-c 7 ^ 


n 


\y,K.,^_„ 2 ) oc 7r(y|K:,’5')7r(/«|^)7r(crf|^_^2) oc Tr{yxt\Kt,ax, Px,crl)TT{al) 

t — 1 X—Xi 

^ ■e^p1-4r^ + jEE(2/ xt {Olx + PxKt))^ J > 


a 


('o.2)np/2+ae-|-l ^ -j ^2 


t X 

The posterior conditional distribution of is thus IG -I- be + ^ EtLi Ex=xi — {<^x + PxKt)) 
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• For crj, we have 


« 7’'(r/|K,^)7r(K|4')7r(cr2|^5r_^2) (xY[Tr{Kt\Kt-i,0,al)Tr(al) 


oc 


(^2)n/2-r5„-ri I 


+ {nt-i + o)y 


The posterior conditional distribution of cr^ is thus IG + ^, bi^ + ^ (^t ~ ■ 


3.4 Forecasting 


The predictive distributions of Un+k^ given y„, are obtained using the MCMC samples as follows. Let L be 
the number of samples remained after burn-in. Then for fc > 1, and for i = 1..., L, we sample recursively 


^n+fe 




Vn+k 


'N 




+ / 3 ' 


^n-\-k 


.0 


2nW 


(15) 


where the samples Kn'^ are obtained from the FFBS procedure. This produces an estimate of '^{yn+klUi-n) = 
/ 7r(y„+j,|K„+fc, ^)7r(K„+fc|K„+fc_i, ’®').. .7r(K„, ^|yi.„) . ..dKn+k and samples from it for fore¬ 

casting. 


4 IMPLICATIONS FOR ANNUITY PRICING 


In this section we aim to quantify the impact of longevity risk on the pricing of annuities, using the mortality 
rates forecasted by the Lee-Carter model in state-space form which is estimated by the Bayesian approach 
described in the previous section. 

4.1 Estimation using Australian mortality data 

The data set consists of Australian female mortality data obtained from the Human Mortality Database 
(http;//www.mortality.org). Since the application is for annuity pricing, we focus on 1-year death rates for age 
60-100fromyear 1975-2011. Figure[T]shows the estimationresults. Here we set = —5,l3xi = 0.2 and as¬ 
sume niQ = 0, Co = 100 (these are the mean and variance of kq used in Kalman filtering), fia = yp = i^e = 0, 
= cr| = CTg = 100, be = doj = 2.1 and be = bui = 0.3. Number of iterations in MCMC is 5000 and the 
burn-in iterations is 1000. We use very vague prior so that estimation is mainly determined by the data and the 
impact from prior is not material. 

4.2 Annuity pricing 


The T year survival probability of a person aged x currently (i.e. f = 0 or year 2012) is determined by 


tPx — lPx+j — 1 — 6 

1=1 1=1 


(16) 


which is a rando m variable since mx+i-^ , for j = 1,..., r, are random quantities forecasted by the 
Lee-Carter model (IDenuit and Dhaenel (l2007h '). Assuming a large enough annuity portfolio, the price of an 
annuity with maturity T year, written for a x-year-old with benefit $1 per year and conditional on the path 


“■IrT 


= ■ ■. ,mx+T-i,T-i), is given by 


T T 

■ B{0, t) E{lT^yx\'mlr) = Y 'r)rPx{mlr) 


T — 1 


r — 1 


(17) 


where 5(0, r) is the r-year bond price, mf. 
survival probability given which is random. 


is th e first r elements of mf.'r , and rPxi'm.i-r) denotes the 
Denuit and Dhaene ( 2007h shows that some bounds of 


xPximf.^) can be computed analytically. BiffisI (2005) evaluates annuity prices allowing for longevity risk 


using financial theory. From an annuity provider’s perspective, what is important, however, is that the annuity 
price is a random quantity depending on the random paths of mf.j,. Moreover, it is important to determine a 
survival curve rPx (as a function of r) in (fTTl i that best captures the mortality experience of the portfolio for risk 
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Survival Curve ( 65 -year-old ) 




Figure 1. (Upper four panels) Posterior mean and 95% confidence interval (Cl) for parameters a., I3\ posterior 
mean and 95% Cl for the latent process k over year 1975-2011; mean and 95% Cl of the predictive distributions 
of log central death rates (?/ 65 j 2 / 70 i 2/75 j 2 / 8 o ) over 40 years forecast. (Lower two panels) Survival curves for 
different ages. 


management purposes. In this regard, we evaluate different quantiles of the annuity price (mf.j,) in Table[T] 
and extract the corresponding survival curves. Note that the forecasted death rate samples are used to produce 
sample paths and hence samples of annuity prices where i = 1,... ,L. The bottom two 

panels of Lig.[T]illustrates the survival curves corresponding to the median, 0.025 quantile and 0.975 quantile 
of the annuity price. Note that a less expensive annuity price indicates smaller survival probabilities. 

Impact of longevity risk. The possibility that the realised survival curve would be different to the survival 
curve assumed for pricing leads to the so-called systematic mortality risk, a.k.a. longevity risk. In Table[T]we 
compare the median, 0.025 quantile and 0.975 quantile of the annuity prices for different ages and maturities. 
We also assume a constant interest rate r = 3% and hence B{0, t) = e"'"’’. Although the price difference 
might appear to be overall small, mis-pricing can be a significant risk when considering a large annuity port¬ 
folio. Lor an annuity portfolio consists of N policies where the benefit per year is B, an under-pricing of 7 % 
of the “correct” annuity price will result in a shortfall of Wi?a^ 7/100 where is the “wrong” annuity price 
being charged with benefit $1 per year. Lor instance, N = 10, 000 policies written to 80-year-old policyhold¬ 
ers with maturity r = 20 years and $20, 000 benefit per year will result in a shortfall of $67 million when the 
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realised survival curve is the one that corresponds to the 0.975 quantile annuity price, while the survival curve 
corresponds to the median annuity price is assumed for pricing (here 7 = 4.1 in Table[T]). Moreover, as shown 
in Table[T] mis-pricing is increasingly more pronounced for older ages as well as for annuity policies having a 
longer maturity. 


Table 1. Annuity price with different age and maturity (T) for female policyholder. Value in bracket ( ) is 
the percentage difference compared to median annuity price. We only consider contracts with maturity so that 
age -I- maturity < 100. 


Maturity (years) 

T = 5 

T = 10 

T = 15 

T = 20 

T = 25 

T = 30 




age 

= 65 



Median 

4.49 

8.18 

11.14 

13.38 

14.88 

15.64 

0.025 Q 

4.48 (-0.2%) 

8.13 (-0.6%) 

11.00 (-1.3%) 

13.10 (-2.1%) 

14.42 (-3.1%) 

15.03 (-3.9%) 

0.975 Q 

4.50 (+0.2%) 

8.22 (+0.6%) 

11.26 (+1.1%) 

13.63 (+1.9%) 

15.31 (+2.9%) 

16.22 (+3.7%) 




age 

= 70 



Median 

4.42 

7.94 

10.57 

12.30 

13.15 

13.41 

0.025 Q 

4.41 (-0.4%) 

7.86 (-1.0%) 

10.37 (-1.9%) 

11.92 (-3.1%) 

12.63 (-4.0%) 

12.82 (-4.4%) 

0.975 Q 

4.44 (+0.4%) 

8.01 (+0.9%) 

10.76 (+1.8%) 

12.66 (+2.9%) 

13.67 (+4.0%) 

14.00 (+4.4%) 




age 

= 75 



Median 

4.31 

7.49 

9.54 

10.52 

10.81 

N.A. 

0.025 Q 

4.29 (-0.7%) 

7.38 (-1.6%) 

9.27 (-2.8%) 

10.12 (-3.8%) 

10.35 (-4.3%) 

N.A. 

0.975 Q 

4.34 (+0.6%) 

7.61 (+1.5%) 

9.80 (+2.8%) 

10.92 (+3.8%) 

11.28 (+4.3%) 

N.A. 




age 

= 80 



Median 

4.08 

6.63 

7.83 

8.18 

n.a. 

N.A. 

0.025 Q 

4.03 (-1.1%) 

6.48 (-2.4%) 

7.57 (-3.4%) 

7.86 (-3.9%) 

n.a. 

N.A. 

0.975 Q 

4.12 (+1.1%) 

6.79 (+2.3%) 

8.10 (+3.4%) 

8.51 (+4.1%) 

n.a. 

N.A. 


5 CONCLUSIONS 

This article explores further the state-space representation of the Lee-Carter model in longevity modeling. We 
derive in details the posterior distributions of the static parameters and the latent process of the model under 
the Bayesian framework via Gibbs sampling. We suggest an identification constraint for the model that is 
particularly suitable for estimation under a MCMC approach. The predictive distributions of death rates are 
used to determine the range of annuity prices. Our results show that the assumption of survival curve has 
significant impact on annuity prices. Annuity written for older age policyholders is particularly vulnerable to 
mis-pricing caused by longevity risk. Extensions of the Lee-Carter model in state-space form and its estimation 
are currently under investigation. 
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